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ABSTRACT 


This  report  is  an  introduction  to  the  Finite  Element 
Adaptive  Research  Solver  (FEARS)  computer  program.  It 
briefly  describes  the  mathematical  problem  solved  by  FEARS. 

It  tells  how  to  prepare  the  input  data  for  FEARS  and  how  to 
interpret  the  output  provided  by  FEARS.  A  composite  shaft 
problem  is  provided  as  an  example. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  David  W.  Taylor  Naval  Ship  R&D  Center's 
(DTNSRDC)  Independent  Research  Program,  Program  Element  61152,  Task  Area 
ZR0140201,  DTNSRDC  Work  Unit  1844-140.  Professor  Ivo  Babuska  of  the  Institute 
for  Physical  Science  and  Technology  (Contract  N001^t)-82M-0743)  and  Mr.  Charles 
K.  Mesztenyi  of  the  Computer  Science  Department,  University  of  Maryland 
(Contract  N00167— 82M-0636)  are  contractors  to  DTNSRDC. 

1.  INTRODUCTION 

The  objective  of  the  Finite  Element  Adaptive  Research  Solver  (FEARS) 
project  wag^to  design  an  innovative  finite  element  solver  and  assess  its 
feasibility  and  computational  efficiency.  Once  completely  developed,  this 
new  finite  element  solver  would  have  the  following  characteristics: 

•  It  would  produce  reliable  information  about  the  total  accuracy  of 
the  finite  element  model.  This  information  would  take  into  account  both  the 
error  due  to  the  finite  element  discretization  and  the  error  of  the  model 
itself  (i.e.,  the  errors  caused  by  the  small  strain  assumption,  the  linear 
constitutive  law,  etc.). 

•  It  would  require  the  user  to  supply  only  the  necessary  physical  infor¬ 
mation,  the  computational  aims,  and  possible  restrictions.  For  example,  the 
user  might  want  to  compute  the  maximum  stresses  or  the  stress  intensity  factor, 
say  with  an  accuracy  to  within  about  5%  and  with  certain  restrictions  pertain¬ 
ing  to  computer  cost,  memory  availability,  etc. 

•  It  would  be  completely  automated,  user-friendly,  and  flexible,  allow¬ 
ing  the  user  to  perform  such  studies  as  sensitivity  analyses  and  evaluations 
of  engineering  data. 

•  It  would  reflect  the  recent  state  of  the  art  in  finite  element  technol¬ 
ogy  as  well  as  current  trends  in  mathematics,  computer  science,  and  engineering. 


The  present  version  of  FEARS  implements  some  of  these  basic  features 
such  as  an  error  estimate  for  the  finite  element  discretization  and  an  adapt¬ 
ive  approach  to  mesh  refinement.  It  al3o  provides  an  option  for  postprocess¬ 
ing.  This  option  will  enable  the  user  to  evaluate  both  the  mathematical  model 
used  and  the  reliability  of  the  engineering  data  and  to  treat  problems  in  frac¬ 
ture  mechanics  among  other  things  when  the  appropriate  postprocessors  are 
implemented . 

The  design  of  the  FEARS  finite  element  solver  is  unique,  and  it  has  capa¬ 
bilities  which  make  it  significantly  different  from  all  present  finite  element 
codes.  Section  2  of  this  report  briefly  discusses  the  main  features  of  FEARS. 
Section  3  describes  the  problem  of  a  cracked  composite  shaft,  formulates  the 
engineering  problem,  and  states  the  problem  in  the  form  directly  solvable  by 
FEARS  (the  F-problem).  Section  4  describes  the  solution  of  the  F-problem  by 
FEARS,  including  preparation  of  the  data  for  the  FEARS  input  of  the  F-problem, 
the  input  in  detail,  and  the  FEARS  output  for  the  F-problem.  Section  5  summa¬ 
rizes  some  results  relevant  to  the  adaptive  aspect  of  FEARS. 

The  basic  mathematical  theory  underlying  the  FEARS  program  was  developed 
at  the  University  of  Maryland.  The  FEARS  computer  program  was  also  designed 
and  implemented  there.  A  numerical  investigation  of  certain  problems  of  spe¬ 
cial  interest  to  the  Navy  is  underway  at  DTNSRDC,  using  the  FEARS  program. 

2.  MAIN  FEATURES  OF  FEARS 

The  current  FEARS  program  has  the  following  features: 

•  It  has  eight  mesh  refinement  procedures  which  are  used  in  an  adaptive 
(or  interactive)  sequence  to  provide  a  solution  with  an  accuracy  specified  by 
the  user.  This  solution  is  the  most  accurate  one  under  the  constraints  (cost, 
time,  etc.)  imposed  by  the  user. 

•  It  provides  reliable  estimates  of  the  error  of  the  finite  element  solu¬ 
tion  with  respect  to  the  (unknown)  exact  solution  of  the  given  boundary  value 
problem. 

•  Its  input  is  both  simple  and  basic.  It  requires  only  the  necessary 
physical  and  geometric  data. 

•  It  has  an  option  for  restarting  the  computation  when  interrupted  by 
the  user  or  by  one  of  the  various  stopping  criteria. 

•  It  provides  a  flexible  frame  for  the  various  postprocessing  procedures. 


FEARS  solves  boundary  value  problems  for  a  linear  self-adjoint  elliptic 
system  of  partial  differential  equations  of  second  order  with  two  field  varia¬ 
bles  in  two  dimensions  when  these  boundary  value  problems  are  formulated  in  a 
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weak  form  on  a  domain  (ICR  .  The  following  discussion  illustrates  the  features 
of  FEARS  in  more  detail . 

a.  DOMAIN  0 

The  domain  0  is  the  union  of  the  interiors  of  a  small  number  (usually 

16)  of  curvilinear  (open)  rectangles  bounded  by  circular  (or  straight)  arcs. 

2 

These  rectangles  are  called  2-D  domains  and  are  denoted  by  D. ,  i  = 

1  i  *• 

The  (open)  bounding  arcs  are  called  1-D  domains  and  are  denoted  by  D.,  j  = 

J 

1,...,N..  The  vertices  of  the  2-D  domains  are  called  0-D  domains  and  are 
denoted  by  D°,  k  =  . 

K  O 

b.  PROBLEMS  SOLVED  (WEAK  FORM) 

FEARS  solves  problems  formulated  in  the  following  weak  form: 


r  r  a  r?£i  j.  vxn  r§£i  .  r^i 

±J  pz\  Ai  LazJ  +  v  Bi  |_3z J  +  UzJ 


Bt  U 
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+  VC, 


U  }  dx  dy 


+  E 
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U  ds  = 


1  |  |  [f]  °(x»y)  +  VTE(x,y)  |  dx  dy 


(2.1) 


+  E  f  e  .  (s)  V  ds 


where 


and  analogously  for  and  V. 

Further 

is  a  4  x  4  symmetric  (constant)  matrix 
is  a  2  x  4  (constant)  matrix 
C.  is  a  2  x  2  symmetric  (constant)  matrix 
Y.  is  a  2  x  2  symmetric  (constant)  matrix 
Di(x,y)  is  a  4  x  1  vector  valued  function 

Ei(x,y)  is  a  2  x  1  vector  valued  function 

Ms)  is  a  1  x  2  vector  valued  function  (s  is  arc  length). 

J 

For  a  discussion  of  the  physical  significance  of  the  matrices,  the 

1* 

reader  is  referred  to  the  User's  Manual.  The  trial  and  test  spaces  used 

allow  us  to  handle  essential  and  natural  boundary  conditions.  The  present 

formulation  allows  us  to  deal  with  materials  which  are  homogeneous  on  every 
2-D  domain  and  with  boundary  conditions  that  are  of  the  same  type  on  every 
1-D  domain. 

c.  ERROR  NORMS 

The  admissible  error  norm  for  the  a-posteriori  estimation  of  the  error 

e  s  u__-u  (where  u_,,  and  u  denote  the  finite  element  and  exact  solution, 
ft  o  FE  O 

respectively)  is  | | | e j | | gp * 


*  The  reference  is  given  on  page  33 


(2.2) 


If  (AE>i  s  Aj,  p  :  1,  Bj  *  0,  Cj  *  0,  and  Yj  =  0,  then  |  |  |.  |  |  |  is  the 
strain  energy  (but  multiplied  by  2).  Because  of  imposed  boundary  conditions 
I  I  I. I | L  is  not  a  seminorm  but  a  norm. 


d.  BOUNDARY  CONDITIONS  AND  RIGHT-HAND  SIDES 

FEARS  provides  a  variety  of  admissible  boundary  conditions  and  right 
hand  sides.  They  can  be  defined  either  in  terms  of  a  class  of  functions 
characterized  by  a  set  of  parameters  or  through  user-provided  subroutines. 


e.  ELEMENTS 

FEARS  uses  elements  of  a  bilinear  type  on  curvilinear  rectangles.  (Tne 

elements  are  bilinear  when  mapped  onto  the  master  square).  Under  certain 

general  conditions  the  asymptotic  accuracy  for  a  smooth  solution  is  of  order 

-1  /2 

h  with  a  rate  of  convergence  of  N  for  a  quasi-uniform  mesh  with  respect 

to  the  energy  norm  (p  =  1).  With  respect  to  the  L  norm  the  asymptotic  accu- 
2  -1 

racy  is  of  order  h  with  a  rate  of  convergence  of  N  (h  denotes  the  size  of 

the  elements  used  and  N  the  number  of  degrees  of  freedom).  The  postprocessing 

allows  us  to  compute  the  functionals  (e.g.,  the  stresses  at  points  of  interest) 

2  -1 

with  an  accuracy  of  order  h  and  a  rate  of  convergence  of  N  (energy  norm). 

The  adaptive  construction  of  the  meshes  leads  to  the  same  rate  of  convergence 
with  respect  to  N  whether  the  solution  is  smooth  or  has  singular  behavior. 

(The  quasi-uniform  mesh  can  result  in  a  rate  of  covergence  N~e  with  e  >  o 
arbitrarily  small  if  the  singularity  of  the  solution  is  very  strong.) 

f.  STORAGE  STRUCTURE 

•  The  storage  structure  of  the  FEARS  program  is  different  from  that  of  the 
usual  finite  'lenient  pro  'ram  due  to  its  adaptive  features.  It  is  based  on  a 
quartic  rooteu  ■  ucture  (father— >4  sons).  Gaussian  elimination  is  used 
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to  solve  the  systems  of  linear  equations  providing  the  finite  element  solution. 
The  process  is  parallel  and  independent  with  respect  to  2-D  domains.  This  fea¬ 
ture  is  reflected  in  the  storage  segmentation  structure.  The  micro-stiffness 
computations,  the  collection  process,  and  the  LU  decomposition  are  performed 
simultaneously  on  the  2-D  domain  level.  The  data  structure  for  the  adaptive 
construction  of  the  meshes  avoids  the  necessity  for  reordering  the  unknowns 
(for  the  purpose  of  elimination)  because  it  is  an  optimal  or  nearly  optimal 
one.  The  unknowns  located  on  the  1-D  and  U-D  domains  are  eliminated  by  a 
sparse  (minimal  degree)  algorithm. 

g.  ERROR  ESTIMATORS  AND  INDICATORS 

The  computation  of  the  error  estimator  is  based  on  the  computation  of 
the  error  indicators  of  each  element,  which  are  computed  independently  on 
every  2-D  domain.  Optimal  meshes  are  characterized  by  the  equilibration 
(equidistribution)  of  the  indicators. 

h.  MESH  REFINEMENT 

FEARS  provides  the  user  with  eight  modes  of  mesh  refinement.  This  allows 
the  user  to  create  arbitrary  meshes  or  have  meshes  adaptively  constructed.  The 
complete  solution  of  the  finite  element  equations  for  the  given  mesh,  including 
computation  of  the  error  indicators  and  estimators,  is  called  the  "long  path" 
solution.  The  "short  path"  solution  provides  the  solution  on  various  sub- 
regions  using  the  nodal  solution  values  outside  of  these  regions  obtained  by 
previous  computation  on  coarser  mesh.  The  "short  path"  solution  on  the  ele¬ 
ment  level  is  inexpesive,  and  the  excessive  use  of  "long  path"  solutions  is  to 
be  avoided . 

i.  DUMP  FILE 

FEARS  provides  various  dump  file  possibilities  which  allow  the  resump¬ 
tion  of  computation  if  processing  is  terminated  by  the  user  or  by  one  of  the 
stopping  criteria. 

j.  OUTPUT  OPTIONS 

FEARS  affords  a  variety  of  output  options  with  respect  to  the  amount  of 
information  to  be  printed.  It  also  provides  information  about  time  used  and 


the  major  characteristics  of  the  adaptive  steps  such  as  statistical  distribu¬ 
tions  of  element  sizes,  sizes  of  the  error  indicators,  etc. 

k.  s  VECTOR 

FEARS  allows  the  user  to  compute  and  print  out  the  five  dimensional 
vector  s  =  (0UTP(1),  0UTP(2) , . . .,0UTP(5) )T 


(2.3) 


and  S  is  a  5  x  6  matrix.  The  values  of  the  vector  s  are  computed  at  the 

T 

centers  of  the  elements.  For  elasticity  problems  s  =  (cs  ,  a  ,  a  ,  u„,  u„) 

1  xx  yy  xy  1  2 


where  a 


aXy  are  the  stresses  and  u ^ ,  u^  the  displacements. 


1 .  POSTPROCESSING 

FEARS  provides  a  frame  for  postprocessing  which  will  allow  the  evaluation 
of  various  functionals  providing  stresses,  stress  intensity  factors,  etc.  with 
higher  accuracy  than  usual  as  well  as  error  estimates. 


3.  PROBLEM  OF  THE  CRACKED  COMPOSITE  SHAFT 


a.  PROBLEM  DESCRIPTION 

The  cracked  composite  shaft  is  shown  in  Figure  3. 1.  The  inner  disk  is 

7 

steel  (Young's  modulus  of  elasticity,  E  =  3.x 10  psi,  Poisson's  ratio,  v  =  .3) 

7 

and  the  outer  ring  is  bronze  (E  =  1.6x10  psi,  v  =  .32).  The  crack  is  at  a  45- 
degree  angle  to  the  horizontal  x-axis  and  extends  halfway  through  the  bronze 
ring  (i.e.,  p  =  r^  +  ,5R).  Plane  strain  behavior  is  assumed.  The  composite 
shaft  is  stress  free  when  the  temperature  of  the  inner  steel  disk  T  =  0  and 
that  of  the  outer  bronze  ring  T  >  0.  Moreover,  on  the  interface  of  the  two 

O 

no  slippage  can  occur.  The  problem  is  then  to  determine  the  stress  state  when, 
after  cooling,  TD  =  Te  =  0.  Obviously  tension  will  occur  in  the  bronze  ring 

D  O 

and  the  crack  will  open  up  inside,  although  both  ends  of  the  crack  stay  closed 
(due  to  non-slippage) .  An  independent  postprocessing  capability  for  FEARS  will 


be  developed  which  can  compute  the  stress  intensity  factors  at  both  ends  of  the 
crack. 


Figure  3. 1  -  Cracked  Composite  Shaft 


b.  F-PROBLEM 

The  problem  must  be  cast  in  a  form  which  FEARS  can  solve,  and  this  FEARS- 
solvable  problem  is  called  an  F-problem. 

The  F-problem:  Find  tne  stress  state  (plane  strain)  of  the  cracked  com¬ 
posite  shaft  when  the  internal  interface  of  the  steel  disk  and  bronze  ring  is 
loaded  with  a  unit  radial  constant  load.  The  accuracy  is  to  be  within  about 
5J  measured  in  the  energy  norm.  (The  energy  norm  of  the  error  is  the  square 
root  of  the  absolute  value  of  the  difference  between  the  energies  of  the  exact 
and  finite  element  solutions  of  the  problem.)  The  F-problem  is  essentially 
equivalent  to  the  cracked  composite  shaft  problem  in  the  sense  that,  when  the 
solution  of  the  F-problem  is  known,  the  solution  of  the  original  problem  is 
easily  obtained.  If  the  outer  bronze  ring  is  allowed  to  cool  (unimpeded  by 
the  inner  steel  disk)  to  the  temperature  T=Q,  the  inner  boundary  is  displaced 


by  an  amount  A  =  Tnr  ,  where  a  is  the  thermal  expansion  coefficient  for 

D  O 

aT  G 

bronze.  On  the  other  hand  a  negative  radial  load  Z  =  - —  placed  on  the 

(l-2v  ) 
s 

interface  creates  a  hydrostatic  pressure  in  the  steel  disk  and  produces  a 
displacement  A  in  the  outer  boundary  of  the  disk.  Therefore,  this  load  in 
the  original  problem  results  in  no  stress  in  the  bronze  ring  and  a  hydrostatic 
stress  state  in  the  steel  disk.  The  additional  stress  caused  by  the  removal 
of  the  auxiliary  load  Z  is  then  computed  by  the  F-problem.  This  removal  is 
effected  by  applying  the  positive  radial  load  Z  to  the  steel-bronze  interface. 
The  combination  of  the  hydrostatic  stress  in  the  steel  disk  and  the  stress 
state  computed  by  FEARS  for  the  F-problem  yields  the  solution  to  the  original 
shrink  fit  problem.  (See  Figure  3.2) 
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Figure  3. 2  -  Scheme  of  the  F-Problem 


4.  FEARS  SOLUTION  OF  THE  F-PROBLEM 
The  FEARS  solution  consists  of  the  following  parts: 

•  Data  preparation 

•  Input 

•  Output 

These  stages  are  common  to  all  FEARS  solutions. 

a.  PREPARATION  OF  THE  DATA 

The  essential  data  of  the  F-problem  consist  of 

•  Geometry 

•  Material  properties 

•  Boundary  conditions  and  loads 

•  Output  computation . 

(1)  Geometry 

The  domain  must  be  partitioned  into  a  set  of  curvilinear  quadrilaterals 
bounded  by  circular  arcs  or  straight  line  segments.  Some  restrictions  are 
placed  on  the  shape  of  the  quadrilaterals  with  respect  to  the  angles  at  the 
vertices  and  the  angles  of  the  circular  arcs.  (See  the  User's  Manual.)  It  is 
essential  for  the  F-problem  of  our  example  that 

•  The  bronze-steel  interface  coincides  with  1-D  domains. 

•  The  sides  of  the  crack  are  distinct  1-D  domains  and  the  ends  of  the 
crack  are  0-D  domains. 

•  Every  2-D  domain  consists  entirely  of  bronze  or  of  steel. 

Figure  4.1  shows  the  selected  partition. 
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Figure  4.1  -  The  Partition  of  the  Domain 
(2)  Material  Properties 

The  bilinear  form  of  Equation  (2.1)  expresses  the  usual  principle  of 
virtual  work.  Accordingly,  in  this  case  B  =  C  =  y  =  0  and 

A  *  E  1-v  0  0 

(1+v) (l-2v) 

0  l-2v  l-2v 

2  2 

n  l-2v  l-2v 

2  2 

v  0  0 

i . 

The  matrices  A  for  steel  and  bronze  are  different.  For  steel  Eg  =  3.  x  10 
psi,  \>s  s  .3,  and  for  bronze  ER  =  1.6  x  10?  p3i  v  _ 
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(3)  Boundary  Conditions  and  Loads 

The  loads  are  defined  through  the  right  hand  side  of  Equation  (2.1). 
Because  the  load  is  on  the  bronze-steel  interface  only,  D  =  E  =  0  and 
the  right  hand  side  consists  only  of  the  term 

E  /^vl  cos  $  v2  s^-n 

3  B] 

where  tne  integration  is  taken  over  all  1-D  domains  on  the  bronze-steel  inter¬ 
face.  <(>  is  the  polar  angle.  (FEARS  allows  direct  input  of  radial  loads.) 

There  is  no  restriction  on  the  u^  except  for  the  specification  of  dis¬ 
placement  for  two  0-D  domains  to  prevent  rigid  body  motion. 

(4)  Output  Computation 

The  desired  stresses  and  displacements  are  given  by  Equation  (2.3)  where 


b .  INPUT 

The  input  for  the  FEARS  program  is  fairly  simple  and  is  divided  into  the 
following  seven  groups: 

•  Job  control  language 

•  Initial  input 

•  Problem  header 

•  Geometry  input 

•  Various  matrices 

•  Initial  subdivision  input 

•  Command  input 
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Table  4. 1  (pages  14-16)  lists  the  input  for  the  cracked  composite  shaft  problem 
discussed  in  Sections  2  and  3  and  displays  these  seven  input  groups  which  are 
described  in  detail  in  the  paragraphs  that  follow.  This  problem  and  its  input 
are  fairly  standard,  but  other  options  are  available  for  different  types  of 
problems.  The  FEARS  User’s  Manual1  should  be  consulted  for  complete  descrip¬ 
tions  of  all  these  options. 

(1)  Job  Control  Language 

The  job  control  language  of  Table  4. 1  is  that  of  the  Cyber  7400  at 
DTNSRDC  under  the  NOS/ BE  operating  system.  The  present  version  of  FEARS 
requires  a  field  length  of  CM23Q000,  which  must  be  specified  on  the  job  card. 
This  field  length  does  not  allow  P4  priority  (prime  express)  runs,  but  one  can 
use  P3  priority  (prime  batch)  if  the  run  will  take  less  than  600  CPU  seconds. 

The  present  version  of  FEARS  resides  in  the  permanent  file  ABSFXYO, 
ID=PUHZ.  The  request  for  the  permanent  file  DUMP  as  well  as  it3  subsequent 
cataloging  under  the  user's  ID  is  needed  because  of  the  DIMP  command  (see 
paragraph  (7)).  FEARS  writes  out  the  computed  results  on  this  file  which 
subsequently  can  be  used  as  input  either  for  restarting  FEARS  or  for  future 
postprocessors . 

(2)  Initial  Input 

The  initial  input  consists  of  two  lines.  The  first  line,  consisting 
of  eight  integers,  specifies  that  the  input  data  and  the  output  produced  by 
each  interactive  step  will  be  printed  on  tape  6.  See  the  User's  Manual1  for 
the  various  options  available  for  the  disposition  of  output.  The  second  line 
consists  of  one  integer.  The  value  zero  indicates  a  start  from  the  beginning. 
For  the  details  of  the  restart  procedure  the  reader  is  referred  to  the  User's 
Manual . 1 

(3)  Problem  Headers 

The  problem  header  consists  of  three  lines.  The  first  line,  consisting 
of  the  two  entries  1,0.,  indicates  that  there  is  no  optional  user-supplied 
subroutine  for  this  problem.  The  second  line,  consisting  of  eight  digits,  is 
the  identification  number  supplied  by  the  user.  The  third  line  consists  of 
one  real  number  p  =  1.0  which  indicates  the  user's  choice  of  norm. 
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S,  1.1. 4,0,  .25 
s2  2. 4. 7,0, .25 
S3  3, 7, 9,0, .25 
s4  4, 9, 14.0. .25 
s5  5.14.17.0. .25 

56  4.17,18.0, .25 

57  7,18.1 .0,. 25 

58  8, 2, 5.0. .5 
s9  9.5,12.0, .5 
S10  10, 12. 19.0, .5 
S17  11,19.2.0.-5 
S1212,3.4.0,0. 

S13  13,6.13,0,0. 

S14 14.13,20.0.0. 
S1515.20.3.0.C. 
S16U,1 .2,0,0. 
S1717,2,3.0.0, 

S18 18,4,5,0,0. 

S19 1 9 , 7 . 8 . 0 . 0 . 

®20 20.8,5,0.0. 
s2i 21.5,6.0.0. 
SM22 ,8.10.0.0. 
S2323, 9. 10.0.0. 
S2424. 10.12.0,0. 
*2525.10.11,0,0. 
®2626, 11 ,12.0,0. 
s27  27,12.13.0,0. 
SM  28, 14, 15. 0.0. 
s29 29, 10, 15, 0,0. 
S30  30,15.16.0,0. 

®31 31,11 ,16.0.0. 
S3232,16.17.0,0. 
Sm33, 16. 19.0.0. 
®34 34. 18. 19, 0.0. 
19.20.0,0. 
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d,  1.2. 1.5. 4 
d2  2. 3. 2. 6. 5 
d3  3,19.13.2,1 
d4  4.20.19,3,2 
d6  5.20.3.13.6 
d6  4. 5. 4, 8, 7 
d7  7,8.7.10,9 
dg  8,5.8.12,10 
dg  9.6.5.13.12 
dio 10 , 10.9,15,14 
d71 11.11 .10.16.15 
<*12 12, 1 2.11.19,16 
d13 13, 13. 12.20.19 
d14 14. 15.14.16,17 
d,g 15, 16.17.19,18 
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(4)  Geometry  Input 

This  group  deals  with  the  geometry  of  the  problem.  The  domain  is  divided 
into  two-dimensional  (2-D) ,  one-dimensional  (1-D),  and  zero-dimensional  (0-D) 
subdomains  as  shown  in  Figure  4.2.  The  0-D,  1-D,  and  2-D  domains  are  numbered 
independently,  starting  with  1. 

(a)  0-D  Domain  (Node)  Data.  The  first  line  contains  the  number  m  of  the 
0-D  domains.  Then  the  lines  n^  i  =  1,...,m  describe  the  l  0-D  domain. 

Each  line  ni  consists  of  six  entries: 

i ,  x ,  y,  b,  u,  v 

i  is  the  index  number  of  this  point,  x  and  y  are  its  coordinates,  b  is  its 
boundary  condition  number  (integer),  and  u  and  v  are  the  prescribed  values  of 
the  enforced  displacement  vector  of  this  point.  The  values  which  b  can  assume 
as  well  as  the  conventions  for  u  and  v,  are  described  in  the  User's  Manual.1 
In  the  present  example  only  points  1  and  4  have  prescribed  displacements. 

(b)  1-D  Domain  (Side)  Data.  The  first  line  contains  the  number  q  of  the  1-D 
domains.  Then  the  lines  s^  i  =  1,...,q  describe  the  ith  1-D  domain  (side). 
Each  line  s^  consists  of  five  entries: 

i,  p,  q,  c,  k 

i  is  the  index  number  of  the  line,  p  and  q  are  the  indices  of  the  end  points 
of  this  1-D  domain,  c  is  its  boundary  condition  number  (integer),  and  k  is  the 
curvature  of  this  1-D  domain.  Conventions  regarding  the  orientation  of  the 
1-D  domain,  the  sign  of  its  curvature,  and  the  values  of  c  are  described  in 
the  User's  Manual.1 

(c)  2-D  Domain  Data.  The  first  line  contains  the  number  j  of  the  2-D  domains 
Then  the  next  j  lines  d^  i  =  1.....J  describe  the  i  2-D  domain. 

Each  line  di  consists  of  five  entries: 

i,  p,  q,  r,  s 
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i  is  the  index  number  of  the  2-D  domain  and  p,  q,  r,  and  s  are  the  indices  of 
its  four  corner  points.  The  ordering  of  the  indices  of  the  corner  points  must 
be  consistent  with  respect  to  clockwise  orientation,  since  the  ordering  for 
the  "master  square"  is 

q'  s' 
p*  r ' 

An  admissible  ordering  is  such  that  by  a  rotation  it  can  be  brought  into 
agreement  with  the  ordering  of  the  "master  square".  For  example, 

r  •  p' 
s'  q' 

is  admissible,  but  neither  of  the  following  two  orderings  is. 

q*  p'  q'  r' 

or 

r*  s'  s'  p' 

See  the  User's  Manual1  for  more  details.  Alternatively  the  FEARS  program 
allows  the  2-D  domains  to  be  input  in  terms  of  their  sides  (which  are  1-D 
domains)  instead  of  their  corner  points.  (The  present  example  specifies  its 
2-D's  by  their  corner  points.) 

(5)  Various  Matrices 

This  group  contains  the  matrices  of  the  bilinear  form  (Equation  2.1),  the 
matrices  of  the  norm  (Equation  2.2),  and  the  output  matrix  (Equation  2.3). 

(a)  Matrices  of  the  2-D  Domains.  FEARS  permits  the  2-D  domains  to  have 
different  sets  of  matrices  which  are  input  in  packages  as  follows: 

The  first  line  of  the  group  consists  of  the  number  of  different  packages 
of  matrices  (£^  =  2).  Then  the  different  packages  are  listed  one  after 
another.  The  first  line  (S,2)  of  a  package  has  the  form 

i .  n  ,  d1 ,  d2, . . .  ,dn 

where  i  is  the  index  of  the  package,  n  is  the  number  of  the  2-D  domains  to 
which  this  package  of  matrices  applies,  and  the  d^'s  are  the  index  numbers  of 
those  2-D's.  Line  consists  of  five  integer  entries  which  indicates  whether 


19 


the  matrices  A,  B,  C,  D,  E  are  present  or  not.  (In  this  example  only  the 
matrix  A  is  present.) 

The  next  four  lines  in  the  package  contain  the  matrix  A  as  given  by  Equat¬ 
ion  (4.1).  This  matrix  A  must  be  entered  according  to  the  numbering  scheme 

1  3  9  11 

2  4  10  12 

5  7  13  15 

6  8  14  16 

The  next  four  lines  describe  the  matrix  A_  for  the  error  norm.  The  reader 

E 

should  also  note  that  in  this  case  A„  =  A.  A.  is  entered  in  the  same  fashion 

h  t 

as  A.  Data  are  read  by  the  FEARS  program  in  "free  format",  i.e.,  without 
reference  to  fixed  FORTRAN  FORMAT  statements.  Accordingly,  A  and  the  other 
matrices  could  have  been  input  using  any  number  of  consecutive  lines  of  data. 

The  next  line  SLh  in  the  package  consists  of  four  numbers  pertaining  to 
the  computation  of  the  norm.  See  the  User’s  Manual  for  the  details.  At  the 
end  of  the  package  is  the  output  matrix  used  to  compute  the  stresses  and  dis¬ 
placement.  The  output  matrix  S,  given  by  Equation  (4.2),  is  a  5  x  6  matrix 
consisting  of  three  stress  components  and  two  displacements.  It  is  input 
according  to  the  following  scheme: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

Again  the  reader  should  note  the  flexibility  afforded  by  "free  format"  as 
regards  matrix  input. 

This  completes  the  input  for  the  first  package  which  describes  the  bronze 
outer  ring  of  the  present  example.  The  input  for  the  remaining  2-D  domains 
pertaining  to  the  inner  steel  disK  is  prepared  in  a  similar  way. 

(b)  Matrices  of  the  1-D  Domains.  This  group  contains  the  input  for  the  line 
integration,  also  broken  up  into  packages.  The  first  line  of  the  group 
consists  of  the  number  of  packages  (£^  =  1).  Then  the  different  packages  are 
listed  one  after  another.  The  first  line  (H^)  of  a  package  is  of  the  form 

i,  n,  Px»  P2 . pn 
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where  i  is  the  index  of  the  package,  n  is  the  number  of  1-D  domains  to  which 
the  package  applies,  and  the  are  the  indices  of  those  1-D  domains.  The 
line  consists  of  two  numbers 

g 

where  g  indicates  whether  line  integration  occurs  with  respect  to  the  bilinear 
form,  and  l  (integer)  indicates  what  type  of  line  integration  is  required  by 
the  right  hand  side  of  Equation  (2.1).  In  this  case  the  polar  load  is  indi¬ 
cated  by  £  i  -1.  The  line  consisting  of  four  entries  relates  the  line  inte¬ 
gration  to  the  bilinear  form.  If  line  indicates  that  no  line  integration 
is  to  take  place,  then  line  8,^  must  contain  appropriate  dummy  values.  The 

line  £  consists  of  four  entries  describing  the  load  vector  with  respect  to 
b 

polar  coordinates  at  the  end  points  of  the  1-D  domain.  For  more  details,  see 
the  User's  Manual. ^ 

(6)  Initial  Subdivision  Input 

This  group,  consisting  of  one  line  of  one  integer  entry,  indicates 
whether  the  initial  finite  element  mesh  has  four  or  sixteen  elements  per  2-D 
domain.  The  present  example  begins  with  the  quartering  of  the  2-D  domains. 

(7)  Command  Input 

The  commands,  which  are  English  words  and  may  be  abbreviated  by  the  first 
four  letters  of  the  word,  are  executed  by  the  FEARS  program  in  the  order  given 
in  the  data  deck.  Most  of  these  commands  must  be  followed  by  one  or  two  lines 
of  "free  format"  entries  which  are  the  arguments  of  that  command.  Only  four 
commands  (ITERATE  (ITER),  PRINT  (PRIN),  DUMP,  and  STOP)  are  considered  here. 

The  words  in  parentheses  are  the  acceptable  shortened  form  of  the  command. 

The  ITERATE  command  directs  the  FEARS  program  to  solve  the  problem  of  the 
data  deck  using  its  automated  iterative  procedure.  The  ITERATE  command  is 
followed  by  a  line  of  four  arguments 

m,  B,  t,  a 

which  specify  the  conditions  under  which  the  program  will  terminate.  (Of 
course,  the  computation  will  be  aborted  if  more  memory  is  called  for  than  is 
available.)  As  soon  as  one  of  the  termination  conditions  is  satisfied,  the 
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iterative  procedure  is  halted  and  the  FEAHS  program  goes  to  the  next  command. 
Here  m  is  the  maximum  number  of  interactive  steps  to  be  performed,  3  is  a 
constant  determining  the  combination  of  long  and  short  paths,  t  (real  number) 
is  the  maximum  allowable  number  of  CPU  seconds  for  the  computation,  and  a  is 
the  required  relative  accuracy  (if  applied  to  positive  definite  problems  with 
respect  to  energy  norms  as  in  the  case  of  elasticity  problems).  0.3  has  been 
found  by  experience  to  be  a  good  value  for  3.  If  the  value  of  6  is  zero,  then 
all  the  iterative  steps  are  long  paths. 

The  DUMP  command  is  followed  by  an  integer  argument  indicating  the  number 
of  the  tape  unit  on  which  this  FEARS  dump  is  to  be  written,  in  this  case  unit 
20.  This  DUMP  command  causes  a  permanent  file  to  be  written  and  enables  the 
user  to  restart  the  computation  from  the  iteration  step  just  before  the  DUMP 
command . 

The  PRINT  command  is  followed  )y  a  line  of  three  integer  arguments 

a ,  b  ,  c 

a  determines  which  type  or  types  of  domains  are  to  be  printed  out,  b  indicates 
what  sort  of  information  is  desired,  and  c  specifies  how  much  information  is 
to  be  printed.  The  most  commonly  used  values  for  these  parameters  are  a  a  -1, 
b  =  -1,  and  c  =  3;  these  provide  all  essential  information. 

The  OUTPUT  command  directs  FEARS  to  calculate  and  print  out  the  stresses. 
It  is  followed  by  a  line  containing  the  integer  1.  Finally,  the  STOP  command 
terminates  the  execution  of  FEAHS. 

c.  OUTPUT 

The  FEARS  program  provides  the  user  with  many  options  for  the  output. 

Only  those  forms  of  the  output  which  the  average  user  would  encounter  are 
described  here.  They  are: 

•  Main  summary  output.  This  is  the  information  provided  by  one 
iterative  step  (either  a  "long  path"  or  "short  path"  solution)  of  the 
automated  solution  procedure. 

•  Total  output.  This  is  the  information  printed  by  the 


PRINT 

-1.-1. 3 


command.  It  contains  all  the  essential  information  concerning  the  mesh  and 
the  solution  at  a  given  step  of  the  automated  solution  procedure. 

•  s  output,  s  is  defined  by  Equation  (2.3).  The  components  of  the 
vector  s  are  denoted  by  OUTP( 1 ),..., 0UTP(5) . 

(1)  Main  Summary  Output 

The  main  output  summary  for  each  iteration  step  begins  with  a  display 
similar  to  the  one  given  here. 
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RESULTING  IN  72  ELEMENTS  WITH  MATRIX  STORAGE 


AND  LONG  PATH  TIME  . 182808E* 


This  display  contains  information  concerning  the  present  mesh  and  number  of 
suggested  elements  for  refinement  in  the  next  iteration  step  along  with  some 
storage  information.  Ordinarily  the  user  will  not  be  concerned  with  this 
display  since  it  may  not  describe  the  actual  refinement  that  is  effected  at 
this  iteration  step.  (This  display  may  not  be  printed  in  future  versions  of 
the  FEARS  program.) 

Ihe  next  display  presents  information  about  the  actual  mesh  refinement 
effected  at  this  iteration  step  as  well  as  information  concerning  the  refine¬ 
ment  of  this  present  mesh  at  the  next  iteration  step.  Information  is  also 
provided  about  the  solution  calculated  for  the  present  mesh.  Here  is  a  sample 
of  such  a  display. 


LONG  PATH  EXECUTION 


*«»a«  FULL  OOMAIN  ***** 

NUMBER  OF  POINTS-  76  MINCER  OF  ELEMENTS-  72 

BORY  POINTS-  55 

ENERGY  NORM-  .601353E-03  ENERGY-  . 3616249787 38E-06 

ERROR  ESTIMATOR-  .  8669490 74 996E-04  RELATIVE  ERROR-  .144167E+00 

MAX. ERROR  INOICATOR-  . 901867667 558E- 0 5  BY  2-0  INDEX-  4  THRESHOLO-  .459925E-P9 

APPROXIMATE  NUMBER  OF  ELEMENTS  TO  eE  SUBOIVIOEO- 
2-0  TO  BE  SU80.  PRESENT  <POINTS» 

10  10  4 

2  2  4  1 

3  0  10  4 

4  2  4  1 

5  0  4  1 

6  0  4  1 

7  0  4  1 

8  0  4  1 

9  0  4  1 

10  0  4  1 

11  0  4  1 

12  0  4  1 

13  0  4  1 

14  0  4  1 

15  0  4  1 

TOTAL  4 

STORAGE  SI2ES= 

2-0  MATRIX,  MAX. CORE  *  112,  TOTAL  =  800,  NO  RECORDS  «  19 

BORY  MATRIX  *  3143 

EXECUTION  TIME* 

SUBOT VISION*  «3730E+00 

2-0  MATRIX  SOLUTION*  .4075E>01 

BORY  MATRIX  SOLUTION*  .2359E*01 

2-0  ERROR  CALCULATION*  . 1585E>01 

BORY  MATRIX  UPDATES  1070  IN  1  RECOROS 


The  number  of  points  is  the  number  of  nodes  in  the  O-D,  1-D,  and  2-D 
domains  which  have  at  least  one  degree  of  freedom.  The  number  of  elements  is 
the  total  number  of  elements  in  the  full  domain.  The  boundary  points  ("BDRY 
POINTS")  entry  is  the  number  of  points  on  the  0-D  and  1-D  domains  which  are 
not  completely  fixed.  The  energy  norm  is  the  square  root  of  the  value  of  the 
energy  when  the  bilinear  form  expressing  the  energy  is  used.  The  error 
estimator  gives  the  error  estimate  for  the  full  domain.  The  relative  error 
is  the  energy  norm  divided  by  the  error  estimator.  The  next  line  gives  the 
maximum  error  indicator,  the  index  of  the  2-D  domain  in  which  it  occurs,  and 
the  threshold  value  for  automatic  refinement.  (All  elements  with  error 
indicators  above  the  threshold  were  refined.)  Then  comes  a  list  of  the  number 
of  elements  in  each  2-D  domain  which  will  be  subdivided  if  the  threshold  value 
is  the  criterion  for  subdivision.  Finally  some  storage  information  and  a 
breakdown  of  the  computing  time  for  that  step  is  provided. 
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(2)  Total  Output 

The  total  output  is  produced  by  the 

PRINT 

-1.-1. 3 

command.  This  output  contains  detailed  information  about  the  0-D,  1-D,  and 
2-D  domains.  Part  of  one  such  display  is  reproduced  here. 


—  3-0  COORDINATES  SOLUTION  ERROR  BOPY  EX 

t  -.400000E+01  0.  0.  0.  0.  0.  3  1 

2  -.20000ae*01  0.  .105294E-07  .677532E-D6  .105294E-07  .6775J2E-08  0  3 

1  0.  0.  19 J643E- 17  . 1 30  282 1- 07  - . 193643E-67  .1332S2E-07  D  0 

4  . 160000E+01  . 3666  06E.01  -. 36970UE - 07  0.  -.369780E-07  0.  2  1 

9  .  800  006E  +  0  0  .1033  03E*  01  -.360329E-07  -.121135E-07  -.  360329E-07  -.121135E-07  0  0 

8  • 800000E+0Q  .80000. E»00  -. 332 789E- 17  .363740E-08  -.332789E-07  .363740E-38  0  0 

7  .260000E  +  01  . 30  39  74f  *01  -.402988E-Q7  . 599245 E-08  -.402988E-07  .599245E-08  0  1 

8  . 20OOOOE+O1  . 2200G0E+01  -.401287E-07  .155474E-08  -.401287E-07  .155474E-08  0  9 

9  . 364575E*  01  «1645  75E*01  -.402101E-07  .174783E-07  -.402101E-07  .174783E-07  0  1 

11  .2870836*01  .8708296*00  -.415106E-37  .161154E-07  -.415106E-07  .161154E-07  0  1 


The  index  numbers  of  the  points  are  given  under  the  "O-D"  heading.  The 
global  coordinates  of  the  points  (x,y)  given  by  the  input  are  shown  under  the 
"COORDINATES"  heading  and  the  computed  solution  values  at  the  point  under  the 
"SOLUTION"  heading.  The  components  of  the  error  of  the  computed  solution  and 
the  exact  solution  (if  the  exact  solution  is  known  and  provided  through  the 
ZMPrRU  subroutine)  are  under  the  "ERROR"  heading.  Otherwise  the  solution 
values  are  printed.  The  boundary  condition  numbers  of  the  points  as  specified 
by  the  geometry  input  are  under  the  "BDRY"  heading.  A  zero  or  a  one  under  the 
"EXP'  heading  shows  whether  that  point  is  internal  (0)  or  external  (1)  to  the 
full  domain . 

Part  of  a  display  of  the  1-D  domain  information  is  reproduced  here. 


---  l-o  I NOE X *  1  B*  0.  E*  1,  1/R«  •  2500Q0E* 00.  FRON  1  TO  4.  PTS*  3.  R-PTS* 

3  POINTS* 

PT  R  LOCAL  SLOBAL  COOP  SOLUTION  ERROR 

1  3  *900000  -.21909E*01  .33466E*01  -.185379E-07  -.931035E-00  -.1S5E-07  -.931E-08 

OERIV  -«67  2 881E-0 1  .105613E*00  0.  0. 

2  3  .290000  -.33900E+01  .21232E*01  -.952604E-08  -.720723E-08  -.953E-08  -. 72 IE- 08 

OERIV  -• 191648E+0  P  ,1180UE*00  0.  0. 

3  3  .750000  -.58997E*00  .39563E*01  -.278484E-07  -.763033E-08  -.278E-07  -.763E-08 

OERIV  -.345002E-01  .149767E+00  0.  0. 
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The  index  number  of  the  1-D  domain  is  clearly  identified.  The  B  value 
is  the  boundary  condition  number  of  the  line  given  by  the  geometry  input.  The 
E  value  shows  whether  the  line  is  internal  (0)  or  external  (1)  to  the  full 
domain.  The  signed  value  of  1/R  is  the  curvature  of  that  1-D  domain  provided 
by  the  geometry  input.  The  indices  of  the  two  endpoints  of  the  1-D  domain  are 
after  "FROM"  (one  index)  and  "TO"  (the  other  index).  (The  orientation  of  the 
1-D  domain  is  indicated  by  the  sign  of  1/R=k.)  The  number  of  points  on  the 
interior  of  that  line  is  the  value  of  "PTS".  The  number  of  regular  points 
on  that  line  is  the  value  of  "R-PTS". 

The  body  of  the  table  provides  information  about  the  points  lying  in 
the  interior  of  that  line.  The  indices  of  these  points  (in  order  of  their 
presentation)  are  listed  under  the  "PT"  Heading.  The  values  of  R  indicate 
whether  the  point  is  regular  (3)  or  not  (0,  1,  or  2) .  A  regular  point  is  a 
nodal  point  which  is  a  corner  point  of  all  the  2-D  domains  on  which  the 
boundary  is  located.  Otherwise  the  point  is  irregular.  (See  the  User's 
Manual. S  The  local  coordinate  of  that  point  on  the  line  is  given  under  the 
"LOCAL"  heading.  The  global  coordinates  (x,y)  of  that  point,  given  under  the 
"GLOBAL  COOR"  heading,  are  the  actual  coordinates  of  the  point.  The  local 
coordinates  are  the  coordinates  of  the  particular  point  on  the  "master"  line. 

The  solution  values  of  the  computed  solution  at  the  point  are  given  under  the 
heading  "SOLUTION".  The  components  of  the  error  between  the  computed  solution 
and  the  exact  solution  (if  the  exact  solution  is  known  and  provided  through 
the  ZPMTRU  subroutine)  are  given  under  the  "ERROR"  heading.  Otherwise  the 
solution  values  are  printed  there. 

The  2-D  domain  information  is  contained  in  a  large  display.  Part  of  the 
information  from  that  display  for  a  sample  2-D  domain  is  displayed  here. 

...  ?-o  INOEX*  1  CORNERS*  2  15  4  NUMBER  OF  POINTS*  31  ELEMENTS*  *6  ENERGY*  .  141  STfc-07 

ELT  SIZES*  0  6  40  0  0  0  0  0,  SMALLER  SIZES*  0,  MAX.  AOJ.  RATIO*  2 

ERROR  INDICATORS,  TOTAL*  .  160 158 E- 09 ,  MAX*  .176372E-10,  MEAN*  .348169E-11,  THRESHOLD*  .1220846-11 

OISTR.SELOH  MEAN  (  35>«  22217667 

0ISTR  .ABOVE  MEAN  (111*  2  9  1  1  1  0  2  4 

OISTR. AMOVE  PREO  (371*  2  0  1  D  0  2  4  26 

TIME,  ASM»OEC*  .4888E+01 ,  8CK*  .5890EA00,  THR*  .1023E*01,  TOTAL*  .6500I.E4L1 

STORAGE,  MEMORY*  360,  AUXILIARY*  560 


The  index  number  of  the  particular  2-D  domain  is  clearly  identified  as  well 

as  the  index  numbers  of  its  corner  points.  The  "NUMBER  OF  POINTS"  is  the 

number  of  regular  points  of  the  mesh  for  this  2-D  domain.  The  value  of 

"ELEMENTS"  is  the  number  of  elements  in  the  mesh  for  this  2-D  domain.  The 

value  of  "ENERGY"  is  the  energy  of  the  solution  on  this  2-D  domain.  The 

values  of  "ELT.  SIZES"  are  the  numbers  of  elements  in  this  2-D  domain  with 

sides  of  length  2_1 ,  2-2 . 2-8,  respectively.  The  value  of  "SMALLER  SIZES" 

—8 

is  the  number  of  elements  in  this  2-D  domain  whose  side  is  smaller  than  2  . 

The  value  of  " MAX. ADJ. RATIO" ,  the  maximum  adjusted  ratio,  is  an  indication  of 
the  maximum  difference  in  sizes  of  two  neighboring  elements  of  this  2-D  domain. 
The  next  line  gives  the  values  of  the  sum,  the  maximum,  the  mean,  and  the 
threshold  for  the  error  indicators  of  the  elements  in  this  2-D  domain.  The 
next  three  lines  present  statistical  information  on  the  distribution  of  the 
indicators.  The  last  two  lines  contain  time  and  storage  information. 

Next  the  nodal  points  of  this  2-D  domain  are  listed;  a  few  lines  of  this 
display  are  reproduced  here. 


PNT 

LOCAL 

COORD 

GL09AL 

COORD 

SOLUTION 

ERROR 

1 

.900000 

.500000 

-•  1 64  317E*  0  1 

•250998E+01 

-. 1425 10E-07 

-.10490  7  E-0  7 

-  .  14  25E-  f  7 

-.lt49E-07 

2 

.Z500C0 

. 250000 

-•  211873E  +  01 

«  1 32  7  01E*0  1 

-.181883E-08 

-,  6550  98E-0  8 

-  .  18 19E- 08 

-•6551L-08 

3 

.290000 

. 750000 

- . 296623E*  0 1 

.  185782E  +  01 

-• 750303E-08 

-.6738 64 E- 08 

-.75  03E-0  8 

-.67J9E-08 

4 

.750000 

.  250  00C 

-.368734E*00 

« 247Z66E+01 

-. 224361E-07 

-.128429E-37 

-.2244F-07 

-.12846-07 

5 

.750000 

.750000 

-.516220000 

, 346172E+01 

- • 259  88  6E-07 

-.867  0  95  E-3  8 

-.25  E9E-07 

-.8671E-08 

6 

.125000 

.125000 

-.21408  0E  *01 

.  692  4326*00 

•452267E-08 

-.179571E-P5 

•4523E-08 

-.1796E-08 

7 

.125000 

.  375000 

-.261654ET31 

« S46306E  +  00 

•  520  266E-09 

-•21264CE-08 

.5203E-09 

-  .2126E-08 

8 

.375000 

.125000 

-•  159887E*0 1 

.1583  086  +  01 

-•449797E-08 

-.  10  7327E-07 

-.44  98E-08 

-.10731-07 

9 

.375000 

. 375000 

-  .1954  17E*  0 1 

. 193487E+01 

-. 824  047E-QS 

-.9060  25E-08 

-, 8240E- 0  8 

-.94606-08 

10 

.625100 

. 125000 

-.8 11 T66E*  0  0 

<  299862E+01 

-•  153457E-07 

-.14583BE-07 

-.15  35E-07 

-.14586-07 

The  index  numbers  of  the 

nodal  points  of  this 

2-D  domain 

are  given 

under 

the  "PNT"  heading.  The  nodal  points  are  listed  in  order  of  creation.  The 
local  coordinates  and  the  global  coordinates  of  this  nodal  point  are  given 
under  the  headings  "LOCAL  COORD"  and  "GLOBAL  COORD",  respectively.  The  com¬ 
ponents  of  the  computed  solution  for  this  nodal  point  are  given  under  the 
heading  "SOLUTION".  The  components  of  the  error  between  the  computed  solution 
and  the  exact  solution  (if  it  is  know  and  provided  through  the  subroutine 
ZMPTRU)  are  given  under  the  heading  "ERROR".  Otherwise  the  solution  values 
are  printed  there.  Only  regular  points  are  listed  here. 
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Finally  the  elements  of  this  2-D  domain  are  shown  in  order  of  decreasing 
error  indicators.  Part  of  one  such  display  of  elements  for  a  given  2-D  domain 


I*.-  4 


is  given  here. 


ELT 

M 

R 

LOCAL 

COORO 

ERROR  INO. 

16 

3 

3 

.125000 

•  875000 

. 176372E-10 

17 

3 

3 

.375000 

.625000 

• 169986E-1Q 

20 

3 

3 

.625000 

.625000 

. 12863 6E-10 

18 

3 

3 

.375000 

.875000 

.  10836  GE-iO 

23 

* 

1 

.875000 

.875000 

.  102  58  5E-10 

21 

3 

3 

.625000 

.875  000 

•  678685E-11 

26 

4 

3 

•  062500 

.062500 

• 553691E-11 

27 

4 

3 

.0  6  25  0  0 

.187500 

• 437586E-11 

59 

4 

3 

.93750  0 

.062500 

•  43139  6E-11 

28 

4 

3 

.187500 

.062500 

• 426943E-11 

29 

4 

3 

.187500 

.187500 

. 379282E-11 

30 

4 

3 

.062500 

•312500 

•  325841E-U 

60 

4 

3 

.937500 

•167500 

. 31940 8E-11 

57 

4 

3 

.812500 

.062500 

. 304  910E-11 

32 

a 

3 

.187500 

.312500 

. 283912E-11 

PREV.ERR. 
• 2548Q1E- 
• 254801E- 
• 266676E- 
• 254801E- 
• 266676E- 
• 266676E- 
• 283420E- 
. 283420E- 
• 277351E- 
.283420E- 
• 263420E- 
• 283420E- 
.277351E- 
•  277351E- 
. 283420E- 


IND. 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 

09 


Ihe  element  indices  are  given  under  the  heading  "ELT".  The  length  of  the  side 

/  u  <  \ 

of  the  element  is  2-v  ”  ,  where  h  is  the  value  of  "H".  The  regularity  condi¬ 

tion  numbers  of  the  elements  of  this  2-D  domain  are  given  under  the  heading 
"R".  A  value  of  3  means  that  the  element  has  four  regular  corner  points; 
values  of  0,  1,  or  2  mean  that  it  does  not.  For  more  details  the  reader  is 
referred  to  the  User's  Manual. ^  The  local  coordinates  of  the  center  of  the 
element  are  given  under  the  heading  "LOCAL  COORD".  The  present  and  previous 
error  indicators  for  this  element  of  this  2-D  domain  are  given  under  the 
headings  "ERROR  IND."  and  "PREV. ERR. IND." ,  respectively. 


(3)  s-Output 

The  OUTPUT  command  prints  out  the  components  of  the  vector  s  given  by 
Equation  (2.3).  Part  of  one  such  display  is  reproduced  as  Figure  4.3. 

The  index  numbers  for  the  2-D  domains  are  given  under  the  heading  "2-D". 
The  index  numbers  for  the  elements  of  a  specific  2-D  domain  are  given  under 
the  heading  "ELT".  If  h  is  the  value  of  "H",  then  the  length  of  the  side  of 
this  element  is  2_v  “  .  The  components  of  the  local  and  global  coordinates 
of  this  element  are  listed  under  the  headings  "LOCAL  COORD"  and  "GLOBAL  COORD" 
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Figure  4.3  -  Components  of  Vector 


The  present  error  indicator  of  the  element  is  given  under  the  heading  "ERR. 

IND.".  It  should  be  clear  that  the  elements  are  listed  in  the  order  of 

increasing  error  indicator.  The  components  of  s  are  defined  by  Equations 

(2.3)  and  (4.2).  The  headings  "0UT(1 )",... ,"0UT(5)"  denote  0  ,  a  ,  a  . 

xx  yy  xy 

u.,u2,  respectively.  These  values  pertain  to  the  centers  of  the  elements. 

5.  SOME  TYPICAL  RESULTS 

Table  5.1  shows  the  dependence  of  the  computed  relative  error  (with 
respect  to  the  energy  norm)  on  the  number  of  elements  used  in  the  long  path 
computations . 


TABLE  5. 1  -  RELATIVE  ERROR  OF  THE  SOLUTION 


No.  of  Elements 
72 
150 
309 
477 


Relative  error 
15.53* 

9.27% 

6.61% 

4.90% 


Experience  shows  that  the  ratio  of  the  true  error  to  the  estimated  error 

never  exceeds  1.2  when  the  error  is  below  10%.  Usually  for  an  error  of  5% 

the  factor  is  1.02-1.06.  (In  theory  the  ratio  of  the  estimator  to  the  true 

error  approaches  1  as  the  error  approaches  zero . ) 

-1/2 

The  table  indicates  a  convergence  rate  of  N  as  it  theoretically 
should  be . 

Table  5.2  3hows  the  number  of  the  elements  in  the  individual  2-0  domains 

TABLE  5.2  -  DISTRIBUTION  OF  THE  ELEMENTS 


Domain 


1 

46 

2 

64 

3 

49 

4 

61 

5 

4 

6 

19 

7 

19 

8 

3M 

Domain 


9 

0 

1 

2 

3 

4 

5 


No.  of  Elements 


The  distribution  of  the  elements  in  the  various  2-D  domains  is  noticeably 
uneven . 

FEARS  provides  a  graphical  display  of  the  meshes  on  the  master  elements 
for  every  2-D  domain.  Figure  5*1  shows  the  meshes  for  2-D  domain  No.  8.  The 
corners  are  numbered  accordingly. 

The  computation  of  the  finite  element  solution  for  the  final  mesh  with 

477  elements  took  about  60S  of  the  total  time.  5X  of  the  final  solution  time 

was  spent  computing  the  error  estimator  and  about  50%  of  this  time  was  spent 

solving  the  linear  algebraic  system  of  equations.  The  time  for  solving  the 

1  65 

system  of  linear  equations  is  roughly  proportional  to  N  *  where  N  is  the 
number  of  degrees  of  freedom. 

The  sequence  of  long  and  short  paths  determined  by  the  ITERATE  command  was 
INIT,  SHORT.  LONG,  SHORT,  LONG,  SHORT,  LONG 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY.  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INIEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


